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Several recent surprises appear dramatically to have improved the likelihood that 
the spin of rapidly rotating, newly formed neutron stars (and, possibly, of old 
stars spun up by accretion) is limited by a nonaxisymmetric instability driven 
by gravitational waves. Except for the earliest part of the spin-down, the axial 
l=m=2 mode (an r-mode) dominates the instability, and the emitted waves may 
be observable by detectors with the sensitivity of LIGO II. A review of these 
hopeful results is followed by a discussion of constraints on the instability set by 
dissipative mechanisms, including viscosity, nonlinear saturation, and energy loss 
to a magnetic field driven by differential rotation. 



1 Introduction 

We review here recent work on a gravitational-wave driven instability that may 
sharply limit the spin of young, rapidly rotating neutron stars. Andersson and 
Kokkotasu have given a comprehensive review of recent work on this r-rnnda, 
instability and its physical implications, following several earlier reviewai3'ESE3. 
The last of these, a general review of the gravitational-radiation driven instability, 
is more pedagogical and not limited to work on r-modes. Included in the present 
review is a survey of recent work on nonlinear saturation and mode coupling not 
included in these earlier articles. 

The first of the surprises mentioned in the abstract was the discovery that the 
r-modes, rotationally restored modes that have axial parity for spherical models, 
are unstable in perfect fluid models with arbitrarily slow rotation. Indicated in 
numerical work by Andersson H, the instability is implied in a nearly. Newtonian 
context by the Newtonian expression fov_thc r-mode frequency E5t!l'L9, and a 
computation by Friedman and Morsink 113 of the canonical energy of initial data 
showed (independent of assumptions about the existence of discrete modes) that 
the instability is a generic feature of axial-parity perturbations of relativistic stars. 
Studies of the viscous and radiative timescales associated with the r-modes 
(Lindblom ct,al.a, Owen et alJE3, Andersson et al. H, Kokkotas and 
Stergioulas E3, Lindblom et al. c3) revealed a second surprising result: The growth 
time of r-modes driven by current-multipole gravitational radiation is significantly 
shorter than had been expected, so short, in fact, that the instability to 
gravitational radiation reaction easily dominates viscous damping in hot, newly 
formed neutron stars (see Fig. || below). As a result, a neutron star that is 
rapidly rotating at birth now appears likely to spin down by radiating most of its 
angular momentum in gravitational waves. (See, however, the caveats below.) 
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Nearly simultaneous with these theoretical surprises was the discovery by 
Marshall et.al. EJ of a fast (16ms) pulsar in a supernova remnant (N157B) in the 
Large Magellanic Cloud. From the pulsar's period and period derivative, and 
from the estimated age of the remnant, the initial period is estimated at less than 
10ms, implying a class of neutron stars that are rapidly rotating at birth. 
Fortifying this conclusion is the belief that accretion-induced collapse of O-Nc-Mg 
and C-0 dwarfs leads to rapidly rotating neutron stars Lj and the likelihood that 
magnetars are rapidly rotating at birth. 

Spurred by these surprises, over fifty authors have worked on aspects of the 
r-modc instability: 



• Examining 
models. lilS 




Studjd; 
emit Ml™ 



rclativistic 



tectability of the gravitational waves they 



• Finding the dominant mechanisms of the effective shear and bulk viscosities, 

mcludirig effects of a suncrflniH int erior and a Solid 

crust, erf 



Finding t he ma x i m r rrrj mode amplitude permitted by nonlinear fluid 
evolution. B@BB@ 



Finding the nonlinear differential rotation associated with r-modes and 
asking whether such^ol^aii dissipates in a magnetic field the energy of an 
unstable r-mode. El 'Hffl'B 



2 The Gravitational- Wave Driven Instability 

All rotating perfect fluid stars are subject to a nonaxisymmctric instability driven 
by gravitational radiation. The instability was found by Chandrasekhar Ej for the 
I = m = 2 polar mode of the uniform-density, uniformly rotating Maclaurin 
spheroids. Although this mode is unstable only for rapidly rotating models, by 
looking at the canonical energy of initiaLdata with arbitrary values of m. 
Friedman and Schutz c2l and Friedman 113 showed that the instability is a generic 
feature of rotating perfect fluid stars, that even slowly rotating perfect-fluid 
models are formally unstable. 

For a normal mode of the form e l ( at + mi p) this nonaxisymmetric instability acts in 
the following manner: In a non-rotating star, gravitational radiation removes 
positive angular momentum from a forward moving mode and negative angular 
momentum from a backward moving mode, thereby damping all time-dependent, 
non-axisymmctric modes. In a star rotating sufficiently fast, however, a backward 
moving mode can be dragged forward as seen by an inertial observer; and it will 
then radiate positive angular momentum. The angular momentum of the mode, 
however, remains negative, because the perturbed star has lower total angular 
momentum than the unperturbed star. As positive angular momentum is removed 
from a mode with negative angular momentum, the angular momentum of the 
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mode becomes increasingly negative, implying that its amplitude increases: The 
mode is driven by gravitational radiation. 

The conclusion, that a mode is unstable if it is prograde relative to infinity and 
retrograde relative to the star is equivalent to requiring that its frequency satisfies 
the condition, 

a{a + mn)<0. (1) 

For the polar f- and p-modes, the frequency is large and approximately real. 
Condition (Q) will be met only if \mil\ is of order |er|, so that for a given angular 
velocity the instability will set in first through modes with large m. 



3 Unstable r-Modes 

Part of the reason that axial modes (r-modes) were not studied extensively for 
neutron stars is that for spherical stars they are stationary, convective currents. ^| 
That is, axial perturbations belonging to an I, m representation of the rotation 
group behave under parity as (— 1) +1 , opposite to Yj m . Because any scalar can be 
written as a sum of Yi m , no scalar perturbation of a spherical star has axial 
parity. Axial perturbations have, to order fl 2 , vanishing Sp,Sp and <5<£>. Only their 
velocity perturbation is nonzero, and it has the form 

Ar l r x VY lm . (2) 

A diagram of the velocity field (due to L. Lindblom) for I = m = 2 is reproduced 
below. 

The fluid elements move along the integral curves of v + <5v, ellipses, to first order 
in the perturbation, as shown in Fig. |^. The restoring force in a rotating frame 
can be regarded as the Coriolis force, leading to a frequency of oscillation 
proportional to angular velocity of the star. 

At the onset of instability, a mode's frequency vanishes. For f-modes and p-modes 
(both polar modes) with low values of I, the frequency of a mode is of order 
cr ~ 1/ ^dynamical ~ w^fi ,J T TTb" TI 5 instability points at similarly large values of the 
star's angular velocityrjO That is, to drag a mode that is retrograde relative to 
the fluid forward relative to infinity, one must have fl ~ c sp h e ricaistarA- I n fact, 
when viscosity is taken into account, polar modes are estimated to be unstable in 
a uniformly rotating star only for il > 0.9Qk, where Q,k is the Kepleri limit on a 
star's rotation, the angular velocity of a satellite at its equator. EJ^oEZIo [] 
The frequency of an r-mode, however is already zero for a spherical star.L_and for 
any nonzero rotation, every r-mode of a perfect-fluid model is unstableU'cB The 
instability in models of slowly rotating, nearly Newtonian stars, follows from the 



"This statement, although formally true, is somewhat misleading. A solution with axial parity to 
the linearized Euler equation is a stationary current, but the nonlinear terms in the Euler equation 
are important once a fluid element has moved a distance of order the radius of the star, and the 
nonlinear solution is time-dependent. 

''This estimate may decrease, perhaps to 0.8 fljc when a general-relativistic computation with 
viscosity is carried out. 
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Flow Pattern for the m = 2 r-mode 
Polar View 




Equatorial View 




Figure 1. The perturbed velocity field of the I = m = 2 r-mode. 



fact that the frequencies El satisfy the criterion (|l|) , 

2(1- l)(l + 2)m 2 n 2 
a(a + mn) = - < 0. (3) 

As noted below, it is likely, but not yet certain that relativistic stars generically 
have a discrete spectrum of r-modes; but a computation by Friedman and 
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Figure 2. Fluid trajectories of the l=m=2 r-mode 



Morsink El of the canonical energy of initial data is independent of the existence 
of normal modes, and it shows that the instability is a generic feature of 
axial-parity fluid perturbations of relativistic stars. 

For stars with general equations of state, the r-modes despdhe the dynamical 
evolution of initial perturbations that have axial paritySO'O In barotropic 
Newtonian models, however, the only purely axial modes allowed are the r-modes 
with I = m and simplest radial behavior .E3 The disappearance of the purely axial 
modes with I ^ m occurs for the following reason.ca Axial perturbations of a 
spherical star are time-independent convective currents with vanishing perturbed 
pressure and density. In spherical barotropic stars, stars for which both star and 
perturbation are governed by a single one-parameter equation of state, the 
gravitational restoring forces that give rise to the g-modes vanish, and they, too, 
become time-independent convective currents with vanishing perturbed pressure 
and density. Thus, the space of zero frequency modes, which generally consists 
only of the axial r-modes, expands for spherical barotropic stars to include the 
polar g-modes. This large degenerate subspace of zero-frequency modes is split by 
rotation to zeroth order in the star's angular velocity, and the corresponding 
modes of rotating barotropic stars arc gencrically hybrids whose spherical limits 
are hybrids of axial and polar perturbations. Because their restoring force is 
rotational (Coriolis), we have referred both to them and to the r-modes as 
rotationaL-modes, and they are called "inertial modes" in the fluid-dynamics 
literature E!l _ . 
To compute these "hybrid" rotational modes, Lockitch and Friedman H expand 
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Figure 3. An r-g "hybrid" rotational mode of a uniform density Newtonian star. The functions 
shown are coefficients of the spherical harmonic expansion (kj) of the perturbed fluid velocity. 



the perturbed 3- velocity, 5v a , in vector spherical harmonics (see also Ref. 

5v a = jr ^WiY l m V a r + ViV a Y l m -iU l e abc V b Y l m V c r\e la \ (4) 

l=rn ^ r ' 

and solve the order fl perturbation equations for the coefficients Ui (r) of the 
axial-parity terms and the coefficients Wi(r) and Vi(r) of the polar-parity terms. 
These coefficients are shown in Fig. || for a particular hybrid mode of a uniform 
density Newtonian model. For a pure r-mode, only one of the coefficients lh(r) 
would be nonzero. 

R-modes of rapidly (and differentially) rotating Newtonian polytropes have been 
recently computed by Karino et al. c3 and references to the earlier Newtonian 
literature are given in the Introduction's bibliography above. 



3. 1 r-Modes of Relativistic Stars 

The r-modeaj of, xotatirig r-rlatjvistic stars were studied fo. 
ntlyffiEZM >BE»BE3 As in the Newtonian case,! 



recc 



first time only 
a spherical 
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Figure 4. Relativistic corrections to the I = m = 2 r-mode of a uniform density star of compactness 
2M/R = 0.316. In Newtonian gravity, only the coefficient U^ir) i s nonzero, i.e. the mode is a 
pure r-mode. However, the relativistic corrections turn the mode into a hybrid. 



barotropic relativistic star has a large degenerate subspace of zero-frequency 
modes consisting of the axial-parity r-modes and the polar-parity g-modes. 
Although barotropic Newtonian stars retain a vestigial set of pure]*-,axjal I = m 
modes, rotating relativistic stars of this type have no pure r-modesJ£3c3 no modes 
whose limit for a spherical star is purely axial. Instead, the Newtonian r-modes 
with I — m > 2 acquire relativistic corrections with both axial and polar parity to 
become discrete rotational modes of the corresponding relativistic models. 
Relativistic mod.es have been computed analytically in a post-Newtonian 
approximations and numerically for slowly rotating polytropes.P°H 66 H 7s l (See Fig. 
IJ.) They have also been studied in numerical time evolutions of rapidly rotating 
relativistic stars using the Cowling approximation (hydro evolution with a frozen 
spacetime metric). L3 Preliminary results suggest that the growth timescale of the 
most unstable mode is largely unaffected by the relativistic corrections, lB differing 
from the post-Newtonian estimatesC3i3 by ^ 10%. 

In the slow-rotation approximation in which they have so far been studied, axial 
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perturbations of non-barotropic stars include, remarkably, a continuous spectrum. 
Kojima E3 shows that the axial modes are described by a single, second-order 
ODE for the modes' radial behavior. (See also the contributions of Kojima and 
Hosonuma to the MG9 r-mode session, APT7). He argues that the continuous 
spectrum is implied by the vanishing of the coefficient of the highest derivative 
term of this equation at some value of the radial coordinate, and Beyer and 
KokkotasB make the claim precise. As the latter authors point out, the continuous 
spectrum they find may be an artifact of the vanishing of the imaginary part of 
the frequency in the slow rotation limit. (Or, more broadly, it may be an artifact 
of the slow rotation approximation.) 

Furthermore, it is possible find discrete normal modes for non-barotropic models. 
This is straightforward in the case of models for which the coefficient in the 
Kojima eqation is nonzero within the fluid. (Non-barotropic uniform density 
models have this property and their r-modes have been computed by a number of 
groups. )BcaE30 On the other hand, it is not so straightforward when the 
coefficient of Kojima's equation vanishes within the fluid - as is the_case for 
certain polytropes with realistic compactness. Ruoff and KokkotasE3 and 
YoshidaEJ have argued that a discrete r-mode simply does not exist for such 
models. However, it would be surprising if a small change in, say, the compactness 
of the star could lead to such a drastic change in the star's physics (the 
disappearance of the r-modes) . It is perhaps more likely that the slow- rotation 
approximation is breaking down in the vicinity of the singular point and that if 
one were to regularize the Kojima equation in this region by including terms-.that 
are higher order in fl, one would be able to find discrete r-mode solutions.aEj 



4 Spin-down and Gravitational Radiation 



Early work suggested that the r-mode instability may sharply limit the spin of 
newly formed, rapidly rotating neutron stars. The radiation may carry off almost 
all of the star's initial rotational energy amoimJbing.±Q few percent of its rest mass, 
and the waves may be detectable by LIGO HcJ-BESO from events out to 20 Mpc. 
The spin-down model on which these tantalizing estimates were based assumed 
that the most unstable r-mode (with multipole indices I = m = 2) would be able 
to grow to an amplitude of order unity before being saturated by some sort of 
nonlinear process, and, as we will see, that assumption has been supported by 
subsequent work. 

Spin-down computations in the context of linear perturbation theory describe a 
competition .between viscosity and gravitational radiation. In these 
calculations ,rHP 8 n3 , c3 the growth time t of an unstable mode (or the damping time 
of a mode stabilized by a viscosity) has the form, 

111 1 

- = + + , (5) 

i *7"GR Tshcarviscosity 'Hjulkviscosity 

with t the e-folding tim£_for each process. (The analysis. sketched below can be 
found in Lindblom ct alcfl; see also Ipser and LindblomEa). When the energy 
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radiated per cycle is small compared to the energy of the mode, the imaginary 
part of the mode frequency is accurately approximated by the expression 

1 1 dE 

t ~ ~2E~dl" 

where E is the energy of the mode as measured in the rotating frame, 



(6) 



2 



( 5 -J +(5$ 



d 3 3 



(7) 



We have, 



dE 



l>2 



\SJ h 



(2r,6a ab Sa* ab + (S668*) , 

where the dissipation due to gravitational radiation @ has coupling constant 

4nG (l + l)(l + 2) 
1 c 2 '+U(Z-l)[(2Z + l)!!] 2 ' 

5a a b and 59 are the coefficients of shear and bulk viscosity; and estimatesBEE 
corresponding coefficients rj and £ are 

2 



T) = 2 x 10 1 



P 



10 15 g-cm 



10*K 
T 



g-cm -s , 



and 



C = 6 x 10 25 



1Hz V 
+ mO J 



P 



10 15 g-cm _ 



(— V 

\UPKJ 



g-cm -s ; 



(8) 

(9) 
I of 

(10) 
(11) 



Polar and axial radiation arise, respectively, from mass and current multipoles 
Dim and J; m , given by the equations, 

Am - / dWpYT 



Im 



Jim — 



dVr l px ■ r x VY, 



(12) 



(/ + l)c. 

The additional factor of v in the current multipoles implies an additional factor of 
v in the radiated energy of axial modes, and hence a smaller expected rate of 
radiation for the same multipole: For a mode of amplitude A = (displacement of 
fluid element)/i?, with R the stellar radius, we have 



dE_ 
~dl 



A 2 M 2 R 2l a 2l+2 



dE_ 
~dt 



A 2 M 2 R 2l+2 a 2l+i . 



(13) 



The extra factor of a 2 in dE/dt for the current multipoles would make polar 
modes dominant if it were not for the fact that their frequencies are small when 
they are unstable. If a newborn neutron star rotates at nearly its maximum 
frequency it is likely initially to be unstable to both polar and axial modes. Most 
of its spin-down, however, should be dominated by the I = m = 2 r-mode. An 
optimistic diagram of this spin-down is shown below in Fig. ^|. (This figure, 
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Temperature (K) 



Figure 5. Critical angular velocity vs. temperature for an n = 1.0 polytrope. Above the solid 
(purple) curve, the star rotates rapidly enough for its fastest growing (I = m = 2) r-mode to be 
unstable, whereas below the curve all modes are damped by viscosity. The dashed curve shows 
the evolution of a rapidly rotating neutron star as it cools and spins down due to the emission 
of gravitational waves. The red curves show the larger critical angular velocities required for 
instability when a crust is present, the upper red curve corresponding to the viscosity when a 
supcrfluid is present. 



aepared by Ben Owen, revises a similar figure in Lindblom, Owen and Morsink 
) The perturbation here is assumed to reach saturation with an amplitude of 
order unity, while still maintaining the character of a linear r-mode. 
The solid (purple) curve shows, for each temperature T, the minimum angular 
velocity above which gravitational radiation can dominate viscosity and drive an 
r-mode instability. Below the solid curve on the right, bulk viscosity due to 
neutrino production in URCA reactions damps the instability. Below the same 
solid curve on the left, shear viscosity damps the instability. The upper red curves 
exhibit the larger shear viscosity from a laminar boundary layer when a crust is 
present (see Sec. 4.1 below). A newborn star that starts with angular velocity Vlx 
follows the dashed trajectory when no crust is present, becoming unstable when 
the temperature drops below about 3 x 1Q 1Q K. The mode quickly reaches its 
saturation amplitude and then radiates the star's angular momentum in 
gravitational waves. Finally, as the star cools and spins down, shear viscosity 
damps the instability. In the diagram, this occurs at a temperature of about 
10 9 _ftT (after about a year, with standard cooling), and at an angular velocity 
below 0.2J7x. As noted in Sect. ^ the amplitude may grow large enough that the 
later part of the nonlinear evolution no longer resembles this linear model. 
This initial scenario is likely to overestimate the duration of the instability and 
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Figure 6. Characteristic gravitational wave amplitude h c (solid curve) compared to the noise 
amplitude h TmB (dashed curves) for LIGO, assuming perfect knowledge of the r-mode waveform. 



the amount of energy radiated, and we discuss below the principal corrections 
that have been considered. 

If the r-modes saturate near unity, a perfect understanding of their waveform 
would allow a high signal to noise ratio for sources beyond 20 Mpc for detectors 
with_the sensitivity of LIGO II, as dramatized by Fig. || (taken from Owen et 
al. Lj with permission of the author). 

In reality, however, the large uncertainty in the waveform substantially reduces 
the prospectsJor detection in the immediate futur£_ A study by Brady and 
T. Creighton O, following the work by Owen et alE^I, sets a lower limit on 
detectability by assuming no knowledge of the source, finding that newly formed 
neutron stars should be detectable by LIGO II with narrow banding out to about 
8 Mpc, with uncertainty allowing a range of perhaps 4-20 Mpc. The Virgo cluster 
would then likely be out of reach, and r-modes could not be detected until more 
sensitive detectors were available. 

A stochastic background of gravitational waves produced by a cosmological 
population of newly formed neutron stars is also likely to-be detected only, by 
interferometers with advanced sensitivities (Owen et al. E3, Ferrari et al. LL3, 
Schneider et al. E9). 

4-1 Role of a Solid Crust 



Beginning with Bildsten and Ushormirsky,0 several authors have considered a 
large ini^eap. in .shear viscosity in a boundary layer near a solid 
crust .□'L3c3u3 , E3'0 Because the melting temperature of the crust is estimated to 
be 10 10 K, with a factor 2 uncertainty, a crust may alter the r-mode instability of 
newly formed stars and will certainly be important in any r-mode instability of 
old stars spun up by accretion. 



If the mode's velocity field vanis 



=s in the crust, it will fall rapidly to zero in a 



boundary layer. In laminar flow,E3 the thickness of the boundary layer (Eckman 
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layer), 

""l/ffi- (14) 

can be estimated by equating the acceleration 8[(dtV-)Wv)] ~ f2fo> of a fluid 
element to the viscous force per unit mass, ~ r]Sv/pd 2 . For a neutron star above 
the superfluid transition temperature, the coefficient of viscosity ?y is of order 
rj = 2 x 10 18 p^5 4 T g -2 gm/cm-s, and d is a few cm. 

As this rough estimate suggests, when a crust is present, dissipation in the viscous 
boundary layer dominates viscous energy loss in Eq. (@) : Because the integrand 
T]Sa ab 5a a b in the layer is enhanced by a factor of order R 2 /d 2 (R the radius of the 
star), the dissipation in the layer is larger than that in the interior by a factor 
R/d. Estimate of the higher angular velocities needed for instability are shown as 
the red curves in Fig. J5. 

Levin and Ushomirskytj point out that this conclusion somewhat overstates the 
case, because the-xrust will participate in an r-mode. Detailed calculations by 
Yoshida and Lee Ej describe the interaction of the stellar r-modes with r-modes of 
the crust, finding a series of avoided crossings. Levin and Ushomirsky find that 
the avoided crossings lead to a wide variation in the fractional drop in Sv across 
the boundary layer, which can range from 1/20 to 1, allowing a possible decrease 
in the boundary layer dissipation by a factor (1/20) 2 = 1/400, compared to the 
estimate above. .— . 
These estimates are based on an assumption of laminar flow, but Wu et alE3 show 
that in a maximally rotating star, a mode whose amplitude exceeds ~ 10 -3 will 
become turbulent. Because the relevant amplitudes are larger than this, theirs is 
the best current computation of viscous damping in the presence of a crust. 
The high boundary-layer viscosity may not significantly alter an r-mode 
instability of newborn neutron stars, because the I = m = 2 r-mode is likely to be 
unstable by the time the temperature has dropped to between 2 x 10 10 A and 
4 x 10 10 A", slightly above the temperature at which the crusLis expected to 
solidify. If this is the case, Lindblom, Owen, and UshomirskyO show that the 
r-mode will prevent a solid crust from forming. Instead, they argue, chunks of ice 
will lie in the outer part of star, with a density adjusted to make the heat 
dissipated by the ice flow balance the mode's gravitational-wave driven rate of 
energy increase. If, as expected, the mode amplitude is above 10~ 4 before a crust 
forms, the ice flow will apparently not significantly affect the r-mode spin-down. 

4-2 Instability in Old Stars Spun-Up by Accretion? 

For old neutron starsLspun up by accretion, the crust is likely to play a crucial 
role. Wagoner's hopeO was that accreting neutron stars would spin up until they 
were unstable (to polar modes, in his paper) and would then radiate angular 
momentum in gravitational w.aves at a rate that balanced the angular momentum 
gained in accretion. Bildsten El and Andersson, Kokkotas, and Stergioulas 13 
considered the mechanism for r-modes in low-mass x-ray binaries (LMXBs), the 
latter suggesting that it might account for the narrow range of observed angular 
velocities. 
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Brown and Ushomirsky E3, however, concluded that low^abserved luminosities 
were inconsistent with the Wagoner mechanism. Levin, Ej reexamined the 
mechanism for r-modes, finding a cycle in which the star is spun up to instability, 
but reaches no equilibrium. Instead, by heating up the star, the mode lowers the 
viscous damping, becomes increasingly unstable, and radiates its energy in a short 
time (less than a year). It is then gradually spun up (for more than 10 6 spars), 
and the cycle repeats. Armed with this revised scenario, Andersson et alo 
suggested that the instability is responsible for the maximum observed angular 
velocity (1.6 ms) of old neutron stars and is consistent with the observed range of 
LMXB spins. 

Despite the large difference between turbulent boundary-layer viscosity and the 
shear_sjiscosity of the neutron-star interior that Andersson et al. □ had used, Wu 
et al.Ea find similar results; that is, like Andersson et al. they find spin-rates 
consistent with those of LMXBs. The maximum angular velocity grows with the 
accretion rate and falls in the narrow range of 490-700 Hz for M /M between 
10~ n and 10~ 8 yr~ x , when the crust is rigid. (If the fractional drop in Sv across 
the boundary layer is only 1/10, the maximum angular velocities fall to 220 - 300 
Hz; and an unexpectedly low melting temperature for the crust will again reduce 
the maximum O.) 

Levin's cycle is qualitatively similar in this analysis, but the spin-down time 
increases to 10 3 yr, and the spin-up time to 10 7 yr. 

In a stuchz_of hypercritical accretion flow onto neutron stars, Yoshida and 
Eriguchi Eil find that an unstable r-mode strongly limits the neutron-star spin. 
The high energy input associated with accretion is large enough that Levin's cycle 
is not seen. These are then a class of systems, possible precursors to compact 
binary systems, in which accreting neutron stars may reach a Wagoner 
equilibrium. _ 
Finally, we should mention a recent preprint by Mendell, Ej examining the 
damping of r-modes by a magnetic field in a viscous boundary layer. The 
damping is large for magnetic fields of order 10 12 G, and, like the boundary-layer 
viscosity, could prevent an r-mode instablity in a newborn star if a crust forms 
before an unstable mode has time to grow. 

5 Nonlinear Calculations 

Much of the very recent work on the r-mode instability has addressed the 
nonlinear evolution of the r-modes. The central issue is whether the instability 
found in idealized models survives the physics that governs a young neutron star. 
We discussed the role of the crust in the last section, and now ask two questions 
related to a mode's nonlinear evolution: Does nonlinear coupling to other modes 
allow an unstable r-mode to grow to unit amplitude? Does the background star 
retain a uniform rotation law as it spins down or does a growing r-mode generate 
significant differential patation? The importance of this last question was 
emphasized by Spruit 113 and by Rezzolla, Lamb and Shapiro E3 who argued that 
differential rotation would wind up a toroidal magnetic field and drain the 
oscillation energy of the r-mode. (See also Rezzolla's contribution to the MG9 
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r-mode workshop APT7). A number of different approaches have since been 
applied to the nonlinear r-mode problem in an attempt to address these questions. 
One notable approach is the direct numerical evolution p£ the nonlinear equations 
describing a self-gravitating fluid. Stergioulas and Font iB have performed 3-D 
general relativistic hydrodynamic_e.volutions in the Cowling approximation, and 
Lindblom, Tohline and Vallisneri ej have performed 3-D Newtonian 
hydrodynamic evolutions with an added driving force representing gravitational 
radiation-reaction, equivalent to that computed previously by Rezzolla et al. lj 
Stergioulas and Font L3 construct an equilibrium model of a rapidly rotating 
relativistic star and add to it an initial perturbation that roughly approximates 
its I = m = 2 r-mode. They then evolve the perturbed star using the nonlinear 
hydrodynamic equations with the spacetime metric held fixed to its equilibrium 
value (the relativistic Cowling approximation). They find no evidence for 
suppression of the mode on a dynamical timescale, even when the mode 
amplitude, A, is initially taken to be of order unity. Because of the approximate 
nature of the initial perturbation, other oscillation modes are excited in the initial 
data. Recall that for a star with a barotropic equation of state, the generic 
rotationally restored mode is not a puce axial-parity r-mode, but an r-g "hybrid" 
mode (Sect. |J). Stergioulas and Font Il3 find that a number of these hybrid modes 
are excited in their initial data with good agreement between the inferred 
frequencies and earlier results from linear perturbation theory H3. In their 
published work, they find no evidence that the dominant mode is leaking its 
oscillation energy to other modes on a dynamical timescale. Instead, a nonlinear 
version of an r-mode appears to persist over the time of the run, about 25 
rotations of the star. In additional runs with amplitudes substantially larger than 
unity however, one no longer sees a coherent r-mode. This may be evidence of 
nonlinear saturation, but further runs adth more accurate initial data will be 
necessary to conclude this definitively o. (See also Stergioulas' contribution to 
the r-mode workshop APT7). 

These conclusions are consistent with preliminary results from-studies of nonlinear 
mode-mode couplings at higher order in perturbation theory 123a. Other r-modes 
of a nonbarotropic star seem to give no indication of a strong coupling to the 
I = m = 2 r-mode unless its amplitude is unphysically large (^4 ^ 30!) E3. Work is 
still in progress on thcjipnlinear coupling of the dominant r-mode to the g- modes 
of nonbarotropic stars and to theJhybrid modes of barotropic stars l£j. 
The results of Stergioulas and Font iB have also been confirmed and significantly 
extended by the calculation of Lindblom, Tohline and Vallisneri. cj In Stergioulas 
and Font's calculation the growth of the unstable r-mode does not occur because 
the spacetime dynamics have been turned off. However, it would be impossible to 
model this growth anyway even in a fully general relativistic hydrodynamic 
evolution, because the timescale on which the mode grows due to the emission of 
gravitational waves far exceeds the dynamical timescale of a rapidly rotating 
neutron star. 

To simulate the growth of the dominant r-mode in a calculation accessible to 
current supercomputers, Lindblom, Tohline and Vallisneri e£I take a different 
approach. They begin by constructing an equilibrium model of a rapidly rotating 
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Newtonian star and add to it a small initial perturbation corresponding to its 
I = m = 2 r-mode. They then evolve the perturbed star by the equations of 
Newtonian hydrodynamics with a post-Newtonian radiation-reaction force that 
drives the current quadrupole associated with the I = m = 2 r-modc. 
By artificially scaling up the strength of the driving force, they are able to shorten 
the growth time of the unstable r-mode by a factor of 4500. In the resulting 
simulation the mode grows exponentially from an amplitude A = 0.1toA = 2.0in 
only about 20 rotations of the star. _ . 
With this magnified radiation-reaction force, Lindblom, Tohline and Vallisncri eB 
are able-to-papfirm the general features of the simplified r-mode spin-down 
models llSEEallJ. In their simulation, the star begins to spin down noticeably 
when the amplitude of the dominant mode is of order unity, and ultimately about 
40% of the star's angular momentum is radiated away. The evolution of the star's 
angular momentum as computed numerically agrees well with the predicted 
angular momentum loss to gravitational radiation. If their model is accurate, 
however, gravitational radiation would not be emitted steadily at a saturation 
amplitude, but would die out after saturation and then reappear as the mode 
regenerates. 

Again, there is no evidence of nonlinear saturation for mode amplitudes A <^ 1. 
The growth of the mode is eventually suppressed at an amplitude A ~ 3.4, and 
the amplitude drops off sharply thereafter. Lindblom, Tohline and Vallisneri 
argue that the mechanism suppressing the mode is the formation of shocks 
associated with the breaking of surface waves on the star. They find no evidence 
of mass-shedding, nor of coupling of the dominant mode to the other r-modes or 
hybrid modes of their Newtonian barotropic model. 

These various studies all provide evidence pointing to the same conclusion: the 
most unstable r-mode appears likely to grow to an amplitude of order unity before 
being suppressed by nonlinear hydrodynamic processes. It is important to 
emphasize, however, that the 3-D numerical simulations have probed nonlinear 
processes occurring only on dynamical timescales and that the actual growth 
timescale for the r-mode instability is longer by a factor of order 10 4 . It is 
possible that the instability may be suppressed by hydrodynamic couplings 
occurring on timescales that are longer than the dynamical timescale but shorter 
than the r-mode growth timescale. Further work clearly needs to be done before 
definitive conclusions can be drawn. Particularly relevant will be the results from 
the ongoing mode-mode coupling studies EaE3. 
Turning to the question of differential rotation, deviations 
law are observed in both of the 3-D numerical simulations 
proposed that differential rotation will be driven by gravitational 
radiation-reactionJ^J as well as being associated with the second order motion of 
the r-mode, itself 13. In a useful toy model, Levin and Ushomirsky lB calculated 
an exact r-mode solution in a thin fluid shell and found both sources of 
differential rotation to be present. 

To address in more detail the issue of whether or not the r-mode ins 
generate significant differential rotation, Friedman, Lockitch and Sa 
calculated the axisymmetric part of the second order r-mode. We work to second 
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order in perturbation theory with the equilibrium solution taken to be either a 
slowly rotating polytrope (with index n = 1) or an arbitrarily rotating uniform 
density star (a Maclaurin spheroid). The first order solution, which appears in the 
source term of the second order equations, is taken to be a pure I = m r-mode 
with amplitude A. 

We find that differential rotation is indeed generated both by gravitational 
radiation-reaction and by the quadratic source terms in Euler's equation; however, 
the latter dominate a post-Newtonian expansion. The functional form of the 
differential rotation is independent of the equation of state - the axisymmetric, 
second order change in v lp being proportional to z 2 (in cylindrical coordinates) for 
both the polytrope and Maclaurin. _ . 
Our result extends that of Rezzolla, Lamb and Shapiro who computed the 
order A 2 differential drift resulting from the linear r-mode velocity field. These 
authors neglect the nonlinear terms in the fluid equations and argue (based on an 
analogy with shallow water waves) that the contribution from the neglected terms 
might be irrelevant. Indeed, for sound waves and shallow water waves, the fluid 
ift computed using the linear velocity field turns out to be exact to second order 
; thus, one may safely ignore the nonlinear terms. However, for the motion of a 
fluid element associated with the r-modes, we find that there is in fact a 
non-negligible contribution from the second-order change in v v . Interestingly, the 
resulting second order differential rotation is stratified on cylinders. It remains to 
be seen whether the coupling of this differential rotation to the star's magnetic 
field does indeed imply suppression of the r-mode instability. 

6 Comments on Future Work 

The overarching questions are whether unstable r-modes limit the spin of either 
young, rapidly rotating neutron stars or of old stars spun up by accretion; and if 
young stars are unstable, how accurately can we characterize the r-mode's wave 
form? 

The nonlinear studies of r-mode instability are still at an early stage, and it is not 
yet certain what limits the amplitude of an r-mode. All computations so far agree 
that there is no saturation until large amplitude is reached, and both numerical 
groups obser-Kft-thc formation of shocks and surface-wave breaking at large 
amplitudes. c3t3 Still unclear, however, is whether sudden damping of the mode 
by wave breaking is an artifact of the numerics or of an amplified radiation 
reaction; the possibility remains that mode-mode coupling enforces a smoother 
limit on the r-mode amplitude. 

It is similarly unclear whether an r-mode will wind up a magnetic field or 
whether, in the presence of a background magnetic field, there is a slightly altered 
mode that does not secularly change the field. 

On the mathematical side, deciding whether the spectrum of rotational modes of 
relativistic stars is continuous or discrete may be a tractable problem and may 
have implications for neutron-star physics. 

More generally, uncertainties in the microphysics (in, e.g., bulk viscosity due to 
hyperon production in the core, or in the coupling of supcrfluid and crust) are 
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large enough that we are unlikely to decide soon whether unstable r-modes play a 
role in the lives of neutron stars. But a discovery of their gravitational waves 
could decide the issue for us. 
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